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We study scaling solutions of the RG flow equation for ^-effective potentials in continuous di- 
mension d > 2. As the dimension is lowered from d = 4 we first observe the appearance of the Ising 
scaling solution and successively the apparence of multi-critical scaling solutions of arbitrary order. 
Approaching d — 2 these multi-critical scaling solutions converge to the unitary minimal models 
found in CFT. 
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Introduction. One of the most valuable properties 
of the functional RG formalism based on the effective 
average action [If is that its use is not restricted to 
any particular dimension, as instead are methods like 
the e-expansion, Montecarlo simulations or conformal 
field theory (CFT) techniques. Often one can obtain 
RG flow equations where d can be considered as an 
external tunable parameter and it is thus possible to 
enquire how properties of the RG "theory space" evolve 
as the dimension is varied. In particular one can observe 
universality classes, i.e. fixed-points of the RG flow, 
emerge or disappear at various threshold or critical 
dimensions. The study of how Z2-universality classes 
depend on d is particularly interesting since it offers the 
possibility of observing how the fixed-points structure 
of theory space interpolates between the upper critical 
dimension d = 4, where there is only the Gaussian 
fixed-point, and d = 2 where one expects, from the CFT 
analysis, an infinite discrete sequence of fixed-points 
describing multi-critical universality classes of arbitrary 
order. Indeed this is what we are able to observe. 

Flow equation for the effective potential. The effective 
average action Tk [ip] is a generalization of the standard 
effective action that depends on the infrared scale k 
The effective potential Vk{f) is found by evaluating the 
effective average action at a constant field configuration 
<p(x) = <p where r^[^] = (J d d x)Vk{f)- In this approxi- 
mation theory space is projected into the infinite dimen- 
sional functional space of effective potentials and the RG 
flow is represented by an exact equation that describes 
the scale dependence of the effective potential 0- In 
terms of dimensionless variables <p — Z k 1 ^ 2 k d / 2 ~ 1 (p and 
V(ip) = k V(0), where Zk is the running wave-function 
renormalization constant^the exact flow equation for the 
effective potential reads 
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with c^ 1 = (47r) d / 2 r(d/2 + I). Note that equation (JTJ 
is a non- linear partial differential equation (PDE). The 
scale dependent anomalous dimension, introduced in ([1]) , 
is defined in terms of the wave-function renormalization 
by fjk — —dt^ogZk] it can be obtained from the dimen- 
sionless effective potential using the following relation Q: 



Vk = C d 
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In ^ ifiQ is the minimum of the dimensionless effective 
potential, i.e. V k ((fo) = 0. 



An interesting feature of equation ([I]) is that it is valid 
for any value of d; in principle it contains information 
about universality classes in any dimension. In particular 
we expect to find non-trivial behaviors in d = 3, where 
equation (Q]) has been extensively studied [J], and in 
d = 2 where one expects a rich behavior, since the CFT 
analysis implies the existence of an infinite number of 
non-trivial RG fixed-points 

Scaling Solutions. Scaling solutions are solutions of 
dtV*(<p) — and correspond to RG fixed-points in the 
functional space of effective potentials. Every scaling so- 
lution, together with its domain of attraction, defines a 
different universality class. 

We start our analysis by considering the case where the 
anomalous dimension is set to zero rj = 0. From ([T]) we 
see that a scaling solution satisfies the following ordinary 
differential equation (ODE): 
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The Z2-symmetry of the effective potential requires that 
its first derivative vanishes at the origin V^(Q) = 0; §5§ 
then implies V*(0) = 1+ c ^ - ) ■ Since equation © is a 

second order non-linear ODE, we need to use numerical 
methods to solve it It's easy to set up the initial 

value problem as a function of the parameter a — V"(0) 
using the two initial conditions just given. 

One immediately observes that for most values of the 
parameter a the solution ends up in a singularity at a 
finite value of the dimensionless field. For every d and a 
we can call this value (p d (a), in this way defining a func- 
tion @. Requiring a scaling solution to be well defined 
for any <p £ R restricts the admissible initial values of 
a to a discrete set {crfj} (labeled by i). One can now 
construct a numerical plot of the function (p d {cr) to find 
the those values where the function (p d (cr) has a 

"spike", since a singularity in (p d {cr) implies that the rel- 
ative scaling solution, obtained by integrating the ODE 
12]) , is a well defined function for every <p £ K. For any d, 
the function <p d (o~) gives us a snapshot of theory space, 
where dimensionless effective potentials are parametrized 
by a and where RG fixed-points, i.e. scaling solutions, 
appear as spikes. By studying (p d {cr) we will be able to 
follow the evolution of universality classes as we vary the 
dimension. 

We can start by studying the function (p d (cr) for d = 4. 
One finds only one spike at c^i = corresponding to the 
Gaussian scaling solution V*((p) = (the singularity at 
a = — 1 is due to the structure of equation ([3]) and does 
not correspond to any scaling solution). The function 
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Figure 1: The function tpf (<r) for (starting from below) d — 4, 3, 2.6, 2.4. In d = 4 one finds only the Gaussian scaling solution, 
represented by a spike at the origin. As the dimension is lowered new spikes emerge from the Gaussian spike, alternatively on 
the left and on the right; these correspond to multi-critical scaling solutions of increasing degree. In particular, in d = 3 one 
observes a spike at negative a, corresponding to the Wilson-Fisher scaling solution; in d = 2.6 the spike at positive a corresponds 
to the tri-critical Ising scaling solution, while in d — 2.4 the newer spike corresponds to the tetra-critical scaling solution. In the 
insets we show the general form of the dimensionless effective potentials obtained by integrating (J3]) using, as initial conditions, 
the positions of the relative spikes, i.e. the values cr*^. 



<Ps( a ) ^ s shown in Figure 1. We find the same qualitative 
result for any d > 4 as expected by the fact that four 
is the upper critical dimension for the Ising (or Wilson- 
Fisher) universality class. 

As we decrease the dimension from d = 4 we observe a 
new spike branching to the left of the Gaussian spike: this 
corresponds to the Ising scaling solution. As we continue 
to lower d the spike moves to the left and for d = 3 the 
function (p\(<j) looks as in Figure 1. As one expects, the 
value of cr* ; 2 at which we observe the Ising spike is nega- 
tive, indicating that the relative scaling solution obtained 
by integrating the fixed-point equation @ is concave at 
the origin. 

For values of the dimension lower than three the func- 
tion <fig(cr) becomes very interesting. Starting at d = 3 a 
new spike branches from the Gaussian one, this time to 
the right: this corresponds to the tri-critical Ising univer- 
sality class, for which this is the upper critical dimension. 
When we reach d — 2.6 this spike is clearly visible, as is 
shown in Figure 1. In d — 2.4 another spike has already 
emerged, this time to the left of the Gaussian spike: this 
corresponds to the tetra-critical Ising universality class. 
In the insets of Figure 1 we show the general form of the 
scaling solutions obtained by integrating equation ([3]) us- 
ing, as initial conditions, the positions of the spikes, i.e. 
the values c*^ for i = 2, 3, 4. One sees that these solu- 
tions indeed correspond to dimensionless effective poten- 
tials with two, three and four minima. 

As we lower further the dimension new spikes emerge 
alternatively on the left or on the right of the Gaussian 
spike. These correspond to multi-critical potentials of 
increasing order. In particular, one observes that the i- 
th multi-critical scaling solution appears at the critical 




Figure 2: The function (ft (a) for values of the dimension 
approaching two; (form bottom) d = 2.1,2.01,2.003,2.001. 
One observes the increasing number of spikes corresponding 
to multi-critical scaling solution of increasing degree. An ar- 
bitrary number of spikes can be observed as d —¥ 2, the only 
constraint being the numerical resolution at which the func- 
tion (ps(cr) is investigated. 



dimension d Cl i] these can be understood as those dimen- 
sions where new operators become relevant as d is low- 
ered. The dimension of the operator tp 2t is 2i(| — 1), 
so the relative coupling A2; scales as d — 2i(| — 1). Set- 
ting this to zero gives d Ct i = -M^; for integer values of 
i = 1,2,3,4,... we find the following sequence of criti- 
cal dimensions d c .; = 00, 4, 3, |, |, |, |, ... Infinity 
is the upper critical dimension for the Gaussian univer- 
sality class, four is the upper critical dimension for the 



Ising universality class, while three is the upper critical 
dimension the tri-critical Ising universality class and so 
on. Note that the critical dimensions d c ^ accumulate at 
d = 2: this is a first hint that in two dimensions there are 
infinitely many different universality classes. By plotting 
the function <Pg(cr) with the appropriate resolution one 
can check that the new spikes emerge precisely at the di- 
mensions d Ci i. It will be intresting to analytically confirm 
these results, as was done for the Polchinski analogue of 
© in 0. 

As shown in Figure 2, when we continue to decrease 
d towards two, the number of spikes of (p\ {a) grows as 
expected. In d = 2 we "lose track" of the multi-critical 
scaling solutions; this is related to the fact that the first 
derivative term on the lhs of §3§ vanishes, making the 
equation exactly integrable. 

Sine- Gordon model. To integrate the fixed-point 
equation ([3]) when d = 2, we recast it in the following 
form: 
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Equation Q can be interpreted as Newton's equation 
where V(<p) 4-» x(t) [lfj; its solution is thus implicitly 
given by 
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bince the "potential" V* - ± log V* in Q is convex and 
bounded from below, the solutions ([3]) are periodic func- 
tions of (p; the period depends on a and the solution with 
zero period is the Gaussian scaling solution. We thus see 
that in d = 2 equation © has a continuum of periodic 
solutions. These are related to the critical sine-Gordon 
model [8j. One can expand the dimensionless effective 
potential in Fourier series V* (</•>) = J2^Li v n,* cos(n/3*y>) 
and check that indeed one finds the fixed-point value 
j8* = V&tt corresponding to the Coleman point. In the 
sine-Gordon model the wave-function renormalization is 
an essential coupling related to the vortex-fugacity; one 
can then use ([T]) and @ to obtain beta functions exibiting 
the Kosterlitz-Thouless-Berezinski phase transition Q . 

Anomalous dimension. In order to follow the evolu- 
tion of scaling solutions down to two dimensions, we turn 
to consider the case of non- vanishing anomalous dimen- 
sion 77 7^ 0. For every scaling solution we will determine 
the allowed values of rj by solving, self-consistently, the 
fixed-point equation 
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0, (6) 



together with © . Since now the second term on the lhs of 
© is non-zero even when d = 2, due to a non- vanishing 
anomalous dimension, we expect to find a discrete set 
of scaling solutions for every d > 2. Following the r\ = 
case, we define the functions (pj.^ 1 (a) and identify the 
discrete values {0^} for which the solutions of ((5]) are 




Figure 3: The anomalous dimensions rji of the first five multi- 
critical scaling solutions as a function of d. 




Figure 4: The anomalous dimensions rji of the first ten multi- 
critical scaling solutions in d = 2. We show the estimates of 
our work (upper curve) together with results of the derivative 
expansion to order d 2 (middle curve) and the exact results 
from the CFT analysis (lower curve). 



well defined for every G K. The functions (p^' v (a), 
for rj ^ 0, turn out to be qualitatively similar to their 
counter-parts (pg(cr) studied previously. 

We proceed as follows: we fix d and we start with an 
ansatz for r); we compute (p^ v (a) from which we find 
the values {ct*'^}; we use them to solve numerically the 
ODE ([6]) to obtain the relative scaling solutions; we esti- 
mate the anomalous dimension by employing ©; we use 
this value as the ansatz for the next iteration until we 
converge to a self-consistent solution of ©. Using this 
procedure we were able to find scaling solutions of ([6]) 
together with the relative anomalous dimensions. In par- 
ticular we started at d = 4 where 77 = and we followed 
the scaling solutions down to d = 2 by making steps 
in d of size 0.1. At every critical d c ,i we were able to 
follow the emerging multi-critical scaling solution. This 
allowed us to calculate the anomalous dimensions rji of 
the multi-critical fixed-points as a function of d. In Fig- 
ure 3 we show the anomalous dimensions of the first five 
multi-critical scaling solutions in the range 2 < d < 3. 
In particular, in d = 2 our fixed-points can be directly 
related to the multi-critical unitary minimal models that 
are found in CFT. 



In d = 2 we compared our estimates for the anoma- 
lous dimensions with those found by Morris [10], using 
the derivative expansion to order d 2 , and with those ob- 
tained from the exact relation r/i — ^ i+1 y^ +i+l , obtained 

by CFT techniques We make the comparison for the 
first ten multi-critical scaling solutions in Figure 4. Due 
to our elementary estimate {2| , our values systematically 
over-estimate both the d 2 and exact values; still it is as- 
tonishing that we correctly reproduce the overall quali- 
tative picture. We stress that we where able to observe 
an arbitrary number of scaling solutions and not only the 
first ten as in [Toj . Our RG analysis predicts an infinite 
sequence of multi-critical scaling solutions of increasing 
degree, with anomalous dimensions decreasing monoton- 
ically with i. It is important to remark that nowhere did 
we assume conformal invariance, but just translation, ro- 
tation, Z2 and RG scale invariance! All our conclusions 
where drawn just from the analysis of the simple ODE 
©■ 

Our analysis can be seen as a RG validation of the cor- 
respondence between Landau-Ginsburg actions and min- 
imal models (Tlj : correspondence that cannot be seen at 
the level of beta functions alone, but requires at least a 
functional truncation of the RG theory space, since only 
in this case it is possible to distinguish spurious fixed- 
points from real ones and thus to disentangle the rich 
fixed-point structure observed near d = 2. This is prob- 
ably one reason why the e-expansion fails to correctly 
detect minimal models in two dimensions (l2^ . 

In view of the results of this paragraph, and of those 
of [hJ, one can expect to reproduce CFT results at 
the quantitative level by employing functional RG 
techniques. Obviously the payback is that the functional 
RG analysis can be done independently and continuously 
in the dimension. 

Critical exponents. Once we have found a scaling so- 
lution, together with its anomalous dimension, we can 
obtain the other critical exponents by studying linear 
perturbations around it. An easy way to do this 0] is 
by expanding the potential in a Taylor series V((p) — 
J2n=i (tfryr^ 2 ™ an d to transform the PDE (gj in a system 
of N tr coupled ODE for the beta functions fan = d t \2n 
of the dimensionless couplings A2„ = k d ~ 2n ^~ 1+ i) \ 2n 
(note that we consider here rj as an input from the pre- 
ceding analysis). RG fixed-points now correspond to so- 
lutions of the following system of algebraic equations: 

ft>„ = n = l,...,N tr . (7) 

As we vary the number N tr of couplings in our truncation 
we usually find more and more solutions to the system 
([7]). Most of these are spurious fixed-points: only with 
the previous analysis of the scaling solutions it is possible 
to understand which of these do correspond to admissi- 
ble solutions of the fixed-point equation (O and which 
instead correspond to singular solutions. Since A2 = cr, 
the values {cr^} that we previously obtained for given 
d and r\ now serve as a "map" to correctly spot the de- 
sired solutions of the algebraic system ([7]). We clearly 
understand that the spurious solutions of |J7J| correspond 
to solutions of ^ which end up in a singularity [fj. 
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10.17 
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Table I: Calculated critical exponents for the Ising universality 
class in two and three dimensions compared to exact results 
Q and "world best" estimates [l3l ]. 



It is useful to proceed as follows. Since the beta func- 
tions 02n are linear in the couplings A2«, we can itera- 
tively solve (UJ to find the couplings A2 n (cr) in terms of 
A2 = cr. Inserting these in 02n+2 = makes it a polyno- 
mial equation in A2 = cr. This can be solved numerically; 
from the solutions we find, we pick the one that best 
approximates the value c*^ relative to the scaling solu- 
tion we are Taylor expanding. Then we obtain the fixed- 
point values for all the dimensionless couplings \2n(&*,i)- 
To determine the critical exponents, we just need to cal- 
culate the stability matrix and evaluate it at \2n(&*,i)'- 
M - dfe„ 

OA 2m „ 

We did this analysis for the Ising universality class. 
In this case the eigenvalues of M nm are such that 
A x (d) < < A 2 (d) < A 3 (d) <■■■ for any 2 < d < 4. By 
varying the number of couplings Ntr m our truncation we 
were able to obtain convergent estimates for the eigen- 
values [lfj. We obtained estimates for the correlation 
length critical exponent V2(d) = — 1/ K\{d) as a function 
of the dimension. We show the results of our analysis 
in Figure 5, together with 772(d) and the other critical 
exponents that we obtained by employing standard 
scaling relations [13j. In two and three dimensions we 
found the results listed in Table 1. Our estimates for 
the exponents ^(d) and 72(d) turned out particularly 
good, both in d — 2 and d = 3. 62(d) and 02(d) are well 
estimated in d = 3 but not so well in d = 2, while the 
value of 02(d) is poor in both cases. The quantitative 
estimate for 772(d), as we observed earlier, is more crude 
due to the elementary ansatz ©. 

Discussion and outlook. In this paper we studied how 
universality classes of scalar theories with Z2-synrmetry 
evolve with the dimension d. Using functional RG tech- 
niques based on the effective average action, we were 
able to follow, continuously with d, the evolution of RG 
fixed-points (represented by scaling solutions) through 
the functional theory space of effective potentials. Even if 
all our analysis was based on the study of a simple ODE, 
we were able to observe a very rich behavior. Above 
four dimensions we found only the Gaussian universal- 
ity class; at d = 3 the Ising or Wilson-Fisher univer- 
sality class appeared; while in two dimensions we found 
a countable infinity of multi-critical fixed-points corre- 
sponding to the unitary minimal models of CFT. More 
importantly, we presented the continuum evolution with 
d of these universality classes; we observed a succession of 
critical dimensions d c ^ (that accumulate at d = 2) where 



2.0 2.5 3.0 3.5 4.0 

Figure 5: Critical exponents for the Ising universality class in the range 2 < d < 4. From top-left 62(d) /5 (empty circles), 72(d) 
(downward triangles), 1*2 (d) (filled circles), 772(d) (squares), /92(d) (upward triangles) and 02(d) (rhombuses). 



new universality classes branched from the Gaussian one. 
The picture that we presented captured all possible criti- 
cal behavior associated to Z2-symmetry in any dimension 
d > 2. We were also able to give quantitative estimates 
for critical exponents as a function of d. 

It will be very interesting to find real physical systems, 
of effective fractal dimension d, where to observe the 
phase transitions related to these universality classes, 
and to compare our critical exponents with experiments. 
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